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ABSTRACT 

This  paper  presents  e  cltss  of  seal-linear  nnaerical 
differentiation  foraulas  which  is  especially  designed  for  functions  with 


steep  gradients.  Basing  on  these,  we  construct  a  seni-lii 


order  difference  scheme  for  solving  the  two-point  singular  perturbation 
problem 

-  e  u* *  +  p(x)  u*  +  q(x)  u  “  f(x),  u(0)  *  u(l)  ■  0 

The  resulting  discrete  semi-linear  system.  AD  *  d  +  0(D),  is  solved  by  a 
simple  iteration:  A  *  d  +  Q(D^  ^),  «  A  Jd  which  is 

2 

convergent  provided  that  the  uniform  mesh  size  h  satisfies  h  £  jjj  e. 

Moreover,  with  the  same  mesh  size  constraint  it  is  shown  that  a  second 

order  semi-linear  scheme  has  one  more  order  of  precision  than  its 

corresponding  central  difference  scheme  for  small  e  and  a  first  order 

error  estimate  l».h  -  ul  <  C  h  is  derived  where  the  constant  C  is 
|  ®  ®  • 

uniformly  bounded  for  all  e>0.  For  the  same  accuracy,  the  present 
nonlinear  scheme  is  much  more  efficient  than  the  usual  linear  centered 
difference  scheme  for  singular  problems.  For  a  linear  and  a  semi-linear 
singular  perturbation  problems,  numerical  results  agreeing  with  the 
above  analysis  are  included. 


HDD 


o 

M  *0  *H 

«K  ©  +J 

2  0  CO 

03  C  V 

i  o  nj  3  'H 

!  H  O  Vi 
n  t-i 

no  s  ij 

I  *-<  M  ft  U 

I  H  fi  -  rs 


n 

®  f . 

I  -c  s 

I  o  O 

o  ^ 

■v  >>  5  « 

c  CO  i-i  . 

O  -H  o 

•h  i — i  -  ;  ® 

P  X>  «  w 

-O  .0  > 

•H  rH  < _ 


!  in  .  > 

pH 

AC  Q  I 


COPY 

IN5.-LCTE0 


2 


ACKNOWLEDGEMENT 

The  author  would  like  to  thank  Profeaaor  Martin  H.  Schultz  for  hie 
helpful  coaaeati  and  Doctors  Dint  Lee  and  Tony  Chan  for  their  help  in 
English  to  improve  this  manuscript. 


I 


3 


1 

L 

j 

* 

1 


I 

I 


I 


1.  Introduction 

Numerical  differential  formulas  play  a  very  important  role  in 
constructing  difference  schemes  of  differential  equations.  Usual 
numerical  differentiation  formulas  based  on  polynomial  approximations 
are  derived  for  smooth  functions  without  large  derivatives,  it  is 
possible  for  these  formulas  to  lead  to  very  poor  results  when  the 
functions  are  not  smooth.  There  are  usually  two  ways  to  avoid  this 
trouble:  refine  the  mesh,  or  use  higher  order  polynomial  interpolation. 
Sometimes  they  are  called  h-version  and  p-version,  respectively. 

The  approach  presented  in  this  paper  is  quite  different.  The  main 
reason  why  the  usual  linear  schemes  lead  to  worse  results  for  problems 
with  large  derivatives,  especially  those  with  singularity,  is  that  the 
usual  numerical  differentiation  formulas  based  on  polynomial 
approximation  is  not  accurate  enough  in  this  case,  e.g.  an  asymptotic 
behavior  near  singularity  is  exponential  .  It  seems  hard  to  get  high 
precision  numerical  differentiation  formula  near  singularity  if  we 
restrict  ourselves  to  use  only  piecewise  polynomial  approximations  or 
other  linear  functional  space.  Hence,  in  this  paper,  we  try  to  look  for 
some  new  numerical  differentiation  formulas  beyond  linear  functional 
space. 

Thus,  two  questions  arise.  The  first  is  how  to  find  such  kind  of 
numerical  differentiation  formulas.  The  second  is  how  to  use  them  to 
construct  difference  scheme  and  how  to  solve  the  resulting  discrete 
non-linear  system  efficiently.  Besides,  we  expect  that  the  resulting 
non-linear  system  should  be  'not  too  far’,  in  some  sense,  from  being 
linear  in  order  to  simplify  both  the  theoretical  analysis  and  practical 
computation.  Therefore,  the  numerical  differentiation  formulas  and  the 
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resulting  schemes,  which  we  introduce  below,  are  called  semi-linear,  or 
in  some  sense  they  are  required  to  he  'weakly*  non-linear. 


The  above  nain  idea  is  based  on  the  author's  previous  work 
([8]-[10]).  The  purpose  of  this  paper  is  to  derive  some  of  these 
seai-linear  numerical  differentiation  foraulas  (section  2)  and  to 
construct  the  corresponding  sani-1 inear  difference  scheaes  for  a  general 
elliptic  singular  perturbation  problem  (section  4):  -eu’ ’+pu'+qu«f , 
u(0)*u(l)«0.  The  same  problem  is  also  used  as  a  test  problem  in  my 
another  paper  [11]  which  extended  the  Ritz-Galerkin  method  from  linear 
subspaces  to  subsets.  A  similar  study  in  two  dimensional  ease  will  be 
discussed  in  [12] .  The  main  idea  of  our  scheme  is  that  we  use  different 
difference  schemes  in  different  subdomains  according  to  the  size  of  the 
discrete  first  divided  difference;  while  it  is  large  the  scheme  will 
choose  a  semi-linear  scheme  automatically,  otherwise  the  scheme  is  the 
same  as  the  usual  linear  scheme,  for  instance,  the  oentered  difference 
scheme.  It  is  shown  that  the  semi-linear  scheme  has  one  more  order  of 
precision  than  the  conventional  centered  difference  scheme  for  the 
singular  perturbation  problem  if  h£2s/llpla,  where  h  is  a  uniform  mesh 
size  in  the  'singular'  subdomain.  In  the  larger  'regular'  subdomain  the 
mesh  size  can  be  used  as  larger  as  we  desire.  Moreover,  for  the  same 
accuracy  the  semi-linear  scheme  costs  less  CFO  time  than  the  linear 
scheme,  being  a  simple  way  to  reduce  the  resulting  semi-linear  system  to 
an  iteration  with  the  corresponding  linear  system.  The  numerical  tests 
presented  in  section  5  match  the  above  analysis  very  well. 

For  convenience,  an  analysis  is  given  in  section  3  for  a  simplified 
model  problem:  -su"+u'*0,  u(0)*0 ,u(l)**l ,  which  was  discussed  by  many 
authors,  e.g.  Ill— [3] . 


2.  ' Sami- linear*  numerical  diffwMttatios  formulas 


Let  u(x)  be  e  function  defined  in  (a,b)  with  a  large  first 
derivative.  Without  loss  of  generality,  let  ns  consider  u(x)  is 
aonotonic  in  the  interval,  and  suppose  that  u  ■  Fx  is  one-to-one  napping 
such  that  x  *  F  1u. 


Denote  x_1  »  a  <  xQ  <  x1  -  b,  and  let  G  be  defined  as  an  indefinite 
integral  of  F,  i.e. 


G(x)  *  /  Fdx,  (2.1) 

By  the  Mean  value  theorem,  there  exist  two  points  z_1  and  z^:  x_^  <  x_4 

<  xQ  <  Xj  <  Xj,  such  that 

[F_1u0,  F"^]  g  «  «<*_!>»  F”lu0]  G  "  (2'2) 


where  [xl,  x2]Y  is  a  notation  of  divided  difference.  Nov  we  look  for  an 
approximate  formula  for  the  first  derivative  at  the  node  x»Xq,  based  on 
the  formulas  (2.2),  as  follows: 

2 

u'(xQ) - UF"1^.  F_1u0]  G  -  [F_1u0.  F~1u_1]  G  ).  (2.3) 

Xl"X-l 

2 

For  the  simplest  case,  take  F  as  the  identity  mapping,  and  G(x)  ■  x  /2, 
then 


®'(*a) 


u(x1)-u(x_1) 


(2.4) 


X1  ”  X-1 

Formula  (2.3)  is  the  same  as  the  usual  central  difference  formula  based 
on  the  quadratic  interpolation.  In  this  case,  it  is  obvious  that  the 
above  nvaerioal  differential  formula  (2.3)  is  an  extension  of  the  usual 
central  difference  formula. 


Now  we  assume  F  to  be  an  admissible  one-to-one  mapping  such  that  G 
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can  be  obtained  from  (2.1)  direetly.  For  any  snob  F,  (2.3)  definee  a 
numerical  formula  for  the  firat  derivative  at  the  node  z  ■  x^. 


Aa  an  example,  let  F  2  f  ■  fr,  where  r  ia  a  real  parameter.  For 
n(x)  >  0  in  (x_1#  Xj),  from  (2.3)  and  (2.1).  the  following  approximation 
formula  ia  obtained 


2r 


(1+rHx^-x^)  u^*  -  u0r 


1+r  1+r  1+r  1+r 

U1  "  u0  Qo  u-i 
{— - - - - - - - )  (2.5) 


r  r 

u0  '  *-l 


where  uQ  -  u(x^),  u^  «  uU^j),  Uj  -  u(x1> . 

Vhen  r  *  j,  zq  "  |  (x^+x_^) ,  the  eorreaponding  formula  to  (2.5)  ia 


,,  .  1.1/2  1/2  w  1/2  1/2  .  1/2  .  ... 
u'(xn)  *  "  tt_i  X  ♦  «0  +  ).  (2.6) 


-1 


Ia  general,  we  have  the  following  error  eatimate: 

k+1  — 1 

Theorem  1*^:  Let  u.  F  a  C  (  x^,  Xj  ),  where  k-3  or  4.  F  u  ia  a 

one-to-one  mapping,  h  *  Xq~x_^  ■  Xj-Xq,  then,  the  remainder  of  the 

nwmerioal  differentiation  formula  (2.3)  equala  to 


i  <  W1^.  F_1.0)  C  - 

h2  d*F_1u  -1 

—  ^f2u"  +  u’2  - T”(5dTS) 

12  **  du2  du 


F  1®_11  G  >  ■  u’Uq)  + 
X)\xmx  +0(hk). 


(2.7) 


Proof ; 

Applying  the  Taylor  expaneion  for  6(y)  upon  x,  one  obtaina 

G(y)  -  G(x)  -  (y-z)G'  (x)  +  (y-z)2  G”(x)/2  +  (y-x)3  G(3)(x)  /3! 
+  (y-x)4  G(4)(x)/4I  +  0((y-z)5) 


Hence 


w<yi*Po*3r-i) 


g(yx)  -  G(y0) 
71  ~  *0 


0(y0)  -  QCy.x) 
yo  -  y-i 


-  y.i 


<V  + 


(yl"y0>  "  (y0_y-l)  (4)  5 

- — —  G  (y0>  +  °((y,-y. > 

4!  v  x  x 


By  aeans  of  rsles  for  finding  derivative  function  in  iaplicit  case,  it 
is  easyly  seen  that 


W(F"1n1,  F_1u0,  F-1u_1) 

h^  d^F~^n  -1 

.  . ,  .  .  d . 2  ,dF  u.-l,  i 

■  to <V  +  -  S'2"  *  *  —2~(  —>  >>, 

12  dn 

+  (Kb**1),  for  k*3 ,4. 


(2.8) 


Q.E.D. 


Corollary  1.  If  the  napping  F  satisfies  the  relation 

d2F-1n  -1 

|_t2n"  +  u'2  - l)l  -  0 

dz  .  2  dn  x*xft 

dn  0 


(2.9) 


then,  the  f omnia  (2.3)  has  an  error  of  third  or  fourth  order  for  k  ■  3 
or  4.  respectively. 


Corollary  2.  If  k  ■  2,  the  error  estiaate  is  of  the  second  order, 
and  the  second  ten  in  the  right  side  of  (2.7)  is  evaluated  at  a  point  ( 

s  (z^.Xj). 

The  d  on  inant  ten  of  the  truncation  error  of  fonula  (2.5)  is 


12  dx  u  x  x0 

This  leads  to  a  fourth  order  nuaerical  differentiation  fonula  if  r  is 


so  chosen  that  the  above  ten  vanishes.  For  instance,  r  can  be  chosen 


2 


« 


8 


? 

it 

i 

/•4 

I 

.  ' 

l 


! 


i 


% 


X 


-1  +  2(  —  - 


-1 


(2.11) 


Both  (2.10)  and  (2.11)  are  implicit  and  they  can  not  be  need  directly 
for  numerical  differentiation.  However,  they  are  neefnl  in  constructing 
a  high  order  difference  scheme  for  differential  equations  ,  even  for 
partial  differential  equations. 


The  limiting  case  of  (2.5)  in  which  r  tends  aero  is  very 
interesting  for  solving  singular  perturbation  problems.  The 
corresponding  formula  with  remainder  term  becomes 


*• (x0)  *  h  { 


u(a1)-u(x0) 


u(x0)-u(x_1) 


) 


2  LogUU^ /u(xQ) )  Log(u(xQ)/u(x_1) ) 

- 4-( 2u"  ~  +  0(h4) 

12  dX  11  X“*0 


(2.12) 


It  follows  from  Theorem  2.1  that  the  formula  (2.5)  is  fourth  order  if 
the  function  u(x)  satisfies  the  non-linear  second  order  differential 
equation 


2uu*  *  -  (l-r)u*2 


bu. 


where  b  “  constant. 


(2.13) 


One  particular  solution  of  (2.13)  is  the  'parabolic*  function 

u(x)  •  (  cx  +  d  )2^*  (2.14) 

where  c,  d  are  constants.  This  means  that  the  'logarithmic  type' 
numerical  differentiation  formula  (2.12).  specifically  designed  for 
solutions  with  large  derivative,  is  also  good  for  smooth  functions,  such 
as  parabolic  function. 

Similarly,  we  can  derive  numerical  differentiation  formulas  for 
higher  derivatives.  For  instance,  there  is  an  analogue  of  Theorem  2.1 
for  second  derivatives. 


Theorem  2.2:  If  u,F  e  C  (x^.Xj),  li  -  xo“x-i  “  xl_x0*  xnd  F  u  is 
one-to-one  napping  in  the  interval,  then 


n'*(x0)  «  -j  UF-1^,  F'1^]  G  +  [F_1n0,  F*"1*^]  G  -  2n()} 
n'  2  d2F_1n  -1 

-i.  dFn,-!,  +0(h2}  (21 

3.1  an  x*xA 

dn  0 


where  G  is  defined  hy  (2.1)  and  u'n  is  given  by  (2.3). 


The  proof  of  Theorem  2.4  follws  the  proof  of  Theoren  2.1,  nsing  the 
Taylor  expansion  and  differentiation  of  iaplicit  functions.  (2.15)  is 


an  extension  of  the  second  order  central  difference  schene,  and  the 
latter  is  only  a  particular  case  of  F  *  I. 
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3.  Model  problem  lul^tii 

In  order  to  give  a  detailed  analysis,  in  this  section  oar  attention 
is  devoted  to  the  following  model  problem  which  is  similar  to  those 
discussed  by  various  authors  cf.  Christie,  Mitchell  (1978),  Zienkiewicz 
and  others  (1978),  Barrett,  Morton  (1980).  Consider 

Lu  *  -su"  +  u'  *  0  in  (0,1),  u(0)  *  0,  u(l)  »  1.  (3.1) 


with  solution  (  Figure  1  ) 


Let  i.  ■  jh  (  j  ■  0,1,..., N;  h  ■  1/N  ),  using  the  conventional 
J 

central  difference  for  the  first  and  second  derivatives  in  (3.1)  leads 
to  the  following  difference  equation 

LjO^  -  -(a+|)U*  2  +  2*0^  +  (|  -  »>uhj+1  ■  0.  1U<N-1.  (3.3) 

with  »  0  and  “  1,  where  a  ■  j  . 


The  exact  solution  of  the  difference  equation  (3.3)  is  given  by 

p<h)J ' 1  _  .tan 

where  p  ■  ^-/2. 


U* 


j 


p(h)N  -  1 


(3.4) 


In  order  to  preserve  the  increasing  monotonic  property  of  the  solution 
(3.2),  it  is  reasonable  to  demand 


a  >  j,  i.e.  h  <  2e.  (3.S) 

Let  E*  -  u(x.)  -11*..  Eh  -  Max  E*  From  (3.2),  (3.4)  and  (3.5),  for 
j  j  J  j 

small  e,  we  have  an  asymptotic  estimate  for  small  h. 


E*  m  xe*/«  -  e1/*  +  e^+x)/«(l-x) 


12s  e2/e  a2 


I. 


in, 


(3.6) 
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with  Eh0  ■  *  0,  and 


0<  E^  eV,"  Eh  -  II*"2  -  7^2- 

Thus,  the  error  will  be  of  zero  order  if  h  is  the  sue  order  as  e. 


(3.7) 


Hence,  in  this  case  the  difference  scheae  can't  converge  for  nail  s, 
unless  the  aesh  size  h  is  excessively  snail.  Meanwhile,  (3.7)  shows 
that  the  aaziana  error  always  occurs  at  the  last  interior  aesh  point. 

The  difficulty  cones  froa  the  fact  that  the  scheae  (3.3)  yields  a  bad 
approxiaate  solution  near  the  layer  boundary  at  the  righthand  end  point. 
As  a  natter  of  fact,  substituting  the  exact  solution  (3.2)  of  (3.1)  into 
(3.3)  yields 


Tt  °j  •  ■  vrJVVi’  *  §  <vrVi> 

2  2 

-  h{  — s  In' '  j  +  J2B<4>j  +  0<l*4)1  +  m*j  +  r*(3>  +  0(h4)) 

-  jj(-su(4)j  +  2u(3,j}  +  0(h5)  -  ~  u(3>j  +  0(h5). 


(3.8) 


where  the  equation  (3.1)  is  used  and 

x/e 


.3  (3).  . 
h  u  (x) 


(e1/8-l).3' 


(3.9) 


Thus,  the  local  truncation  error  near  x«l  is  only  of  zero  order  if  s  is 
the  sue  order  of  h. 


Therefore,  it  is  why  we  need  to  look  for  a  schue  with  better 
approximation  near  the  layer  boundary.  To  construct  such  a  scheae  we 
try  to  use  the  results  described  in  last  section. 


Froa  (3.2),  we  observe  that 


u(3)(x) 


_  d_  u' 

(e1/e-l)«3'  dx  n 


-  1 
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clearly  both  terns  of  the  above  formulas  are  very  large  vith  the  sane 
sign  near  the  layer  boundary,  but  their  difference  is  very  snail  for 
snail  t  if  x  is  close  to  the  right  end.  in  fact 


(3).  , 
u  (z) 


d  u'2  .  -3  -1/e 

- - >  e  e 

dz  u 


as  z  — >  1. 


Hence,  if  we  could  find  a  scheme  of  (3.1)  such  that  the  doninant 
(3) 

coefficient  u  in  the  truncation  error  formula  (3.8)  could  be  instead 
2 

by  u^(z)  -  4-  * — *  it  would  greatly  reduce  the  approzimation  error 
dz  u 

near  the  layer  boundary.  Therefore,  it  reminds  us  to  adopt  the 
'logarithmic  type'  numerical  differential  formula  (2.12)  near  layer 
boundary. 


Thus,  the  difference  scheme  (3.3)  is  changed  now  as  follows: 

-  -(a  +  |)Dhj_1  +  2au\  +  (|  “  for  2j  <  N  +  1 

L/j  -  -'•4,0Vi  *  "Y  "V- 

for  N+l  <  2 j  .<  2N  (3.10) 


where 


Q(V1*Dj,V1) "  • 

0(zj)-0(zj_1) 

Log(0(zj)/D(zj_1)) 


-  ( 


Log(U(zj+1)/0(zj)) 


with  D 


0  and  U1 


N 


-  » 
1,  *  ■  £• 


(3.11) 


The  looal  truncation  error  of  the  scheme  becomes 


if  2j  <  N  +  1 


Lh  *  ■  ( 


(3.12) 
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i!  /  <3>_  d_  u •*. 

12  {  dx  u  }  j 


if  Nfli2ji2N 


<e1/e-l)e3 


12h“3Trh . 

J 


if  2j  <  N  +  1 


(3.13) 


.  x/e  , .23 
( e  -1 )  e 


if  N4-li2ji2N 


where  x  ■  jh. 


Denote  the  two  right  hand  sides  of  (3.13)  by  Rl(x)  and  R2(x), 
respectively.  Let 

»'*>  *  SSfr  - 177— • 

©  "A 


Sinoe 

#l/(2a)-:l 

D*(x» 0.  D(0)-0,  D(l)-e1/#-l,  D(jb  -  --77--  .  <  1. 

1  el/ize)+1 

the  solution  of  the  equation  D(x)  ■  1  is  equal  to 

Then  the  ratio  function  Bl(x)/B2(x)  increases  frca  xero  to  a  very  large 
nuaiber  in  the  interval  [0,  1],  and  the  asyaptotic  solution  of  Rl(x)  “ 
S2(x)  is  equal  to  ^  +  o(e),  for  saall  s.  This  is  reasonable  to  choose 
as  the  separate  point  of  the  two  different  scheae  in  (3.10).  and  now 
the  aaxiaua  local  truncated  error  of  (3.10)  is  given  by 

3  2  e1/(2e)  3 

Trk  ■  Max  Trh  -  —82(1/2)  -  - - -  h_  -l/ (2e )  -3 

mx  ir  j  12  j  1/(2,)^  2  126  * 


It  is  worth  to  note  that  there  is  a  striking  contrast  between  the 
truncation  errors  for  the  revised  soheae  (3.10)  and  the  original  linear 


*J|M 


L5 


I 

I 
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scheme  (3.3).  In  the  linear  scheme  (3.3)  the  maximum  truncation  error 

has  no  power  of  the  mesh  size  h.  Bnt  the  truncation  error  of  the  revised 
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scheme  (3.10)  has  an  h  term  with  a  coefficient  which  tends  to  0  as  e 
does.  Naturally,  we  expect  the  revised  scheme  (3.10).  which  we  will  call 
a  semi-linear  scheme  based  on  the  numerical  differentiation  formula  in 
section  2,  to  give  much  better  numerical  results  for  the  Model  equation 
(3.1)  than  (3.3)  does. 


However,  the  first  important  question  arisen  is  how  to  solve  the 
semi-linear  system  (3.10).  Denote  the  resulting  system  by 


A  U  -  d  +  Q(  D  ) 


(3.15) 


where  A  and  d  coincide  with  the  corresponding  linear  scheme,  and  the 


non-linear  term  Q(U)  arises  from  using  the  semi-linear  scheme  for  x  > 


Let 


2a  -(a-1/2) 

-(a+1/2)  2a  -(a-1/2) 


-(a+1/2)  2a  -(a-1/2) 
-(a+1/2)  2a 


n-1 


D  *  Det(A  ),  fi  «  1 1 .  then 

n  n  n  ^ 


2  1, 


Dn  “  2‘Dn-l  "  U  i)Dn-2'  V1'  D1  “  2*' 


.  P  ,  *  0. 


2  '  r-l 

2a-(a  -l/4)p  . 

n— a 


h  *  *  5i$f  -  W- 


»  ^n  i  M  11*  L  ■  * 

n  7+m- 


The 


.-1 


inverse  matrix  A  can  be  easily  found: 


(3.17) 
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-1  ,  _1  * 

A  *  (a  ), 

*»  J 

where  a"1.  -  {  (a+l/2)i_jD.  ,D  ./D  if  i>j 

i»j  j-1  n-i  n 


(a-l/2)j_iD.  ,D  ,/D  if  i<j 
i-1  n-j  n 


(3.18) 


We  seed  to  introduce  the  following  Leases: 


1  .-1 


For  *  2  j,  A  is  e  positive  aetrix  given  by  (3.18). 


Moverover 


1  ‘"S.j-i  u  111 


i .  ^  if  1  <  J 


I*'1!,  i  (K-DVl  < 


Suppose  that  a,b,  and  e  are  positive,  let 


(3.19) 


r(s,b)  ■ 


a+b  b  -  a 
2  “  Log(b/a) ’ 


Q(a,b,c)  a  r(b,e)  -  r(a,b) 


(3.20) 


1  (b  -  a)/b 

Qa  *  at  m  2 - + - 2  ’ 

Log(b/a)  (Log(b/a)) 

dQ  i  1  <«  ”  b>/c 

Q  |  ZJZ  m  =  -  —  +  _____ 

C  *C  *  Log(c/b)  (Log(e/b)}^ 


(3.21) 


s-S? 


b  -  a 


(b  -  s)/a  1 


(c  -  b)/b 


Log(b/a)  (Log(b/a))  Log(c/b)  (Log(c/b))‘ 


The  tens  Log(b/«  j  is  a  Generalised  Mean  defined  by  Jiachang  Sun  [9] 
between  the  arithaetic  Means  and  Geoaetric  Means  of  a  and  b: 


L*MM  3*2:  If  a,b  >  0,  then 


(ab)1'2  <  a  -  b  ,  *+b 
Eogla/b)  <  ~T 


(3.22) 


where  the  equality  holds  if  end  only  if  a*b. 


In  feet,  integrating  t  froai  0  to  1  on  the  both  sides  of  the  following 
inequalities:  a1  *b*  i  (l+t)a  +  tb,  it  leads  to  the  righthand  of 
(3.22),  and  the  lefthand  can  be  obtained  by 


/(j1  (J)4"172  dt  «  /0172  {(J)*"172  +  (j)t_172}  dt  2  1. 


Corollary.  For  a,b  >  0 


0  .<  r(a,b) 


a+b 

2 


»  -  b  <  1 
Log(a/b)  2 


<b1/2-.1/2)2. 


(3.23: 


There  exist  a  series  of  inequalities  for  the  function  Q  and  its 
derivatives,  based  on  (3.22). 


£ea»  If  a,b,  and  c  are  positive,  then 


...  1  ,.1/2  1/2.2  .  A/  .  ,  .  1  ,  1/2  .1/2,2 

(i)  -  j  (b  -a  )  4  Q(a,b, c)  i  j  (c  -b  )  . 


(ii)  If  o2b2*>0  and  b  iac,  then  Q(a,b,e)  2  0  with  iff  c  *  a, 


(iii)  0  i  |  ((J)174  "  ll  i  Jf  *  \  -  1).  i*  b2a>0 

-  \  (1-(J)172)  i  (l-(J)l72))i  0,  if  a2b>0. 

0  (1-(|)172)  1  |2  i  |  (l-£),  if  c2b>0 

(J-l)  iff  i  ((J)1/4-l)i  0,  if  b2c>0. 

Now  we  look  for  an  estiaate  of  ||A~2J(Q(C) )M.  where  J(Q)  is  the  Jacobi 
aatrix  of  Q.  First,  we  prove  a  very  interesting  relationship  for  the 
scheae  (3.10) 


(J(Q(u))  u)j  -  (Q(u))j,  if  j  <  N-l.  (3.24] 


Usually,  (3.24)  is  only  true  for  Q(a)  *  Au,  where  A  is  a  aatrix,  it 
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shows  again  that  the  seai-linear  ten  Q(u)  is  really  not  far  froai 
linear.  In  fact,  froai  (3.11)  and  (3.21) ,  (J(Q(u))  u)j  *  0  if  j  <  n. 
where  n  “  [N/2].  And  for  n  <  j  <  N-l 


»),  -  5^ Vl  +  ISpj  * 


»u  ,>-•(.  > 


_  {  _ il± _ i _ 1 _ - } 

Log(n(Xj+1) /n(»j ) )  Log (n(x^ ) /u(Xj_j) ) 


Q.(u) 

J 


Hence 


J(Q(n))  n 


(0i...#0f Q  »  . 

n 


,0N-2'Q  N-l}' 


where 


Qj  *  rj+l/2  "  rj-l/2* 


n.  -  #, 


v  m  -a  -e  -  *  . . . a-  ■  . 

j”1/2  2  Log(u./u.  ) 

*  ®Vl 

«k-i- Vi-  V  ( 

Using  the  shore  results,  we  derive  an  upper  bound  of  a  non  of  the 
■atrix  A  *J(Q(u)).  A  straightforward  computation  yields 


(3.25) 


*1  *  IA'1J«K»»»>J  ’  -  Vn-!  ‘ 


jLl‘  *  ‘.J  '  “  i.J-1  lrJ-l« 


”  a  i,nrn-l/2 


(3.2 6) 


for  n  i  i  i  N-l,  where 


*N.N-  1  "  rN-l/2 


cN-l/2 


©i  »  u  i 


—  -  -i 

,N-1  ®N,N-1  +  j£i+1{  °  i.j- 


1  ‘  °  1i,j  )rj-l/2 


a  i#j  ”  a  i.j-1  ,rj-l/2  “  a  i,nrn-l/2* 

'#♦1 
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Suppose  >  0,  ■  1,  using  (3.23)  end  (3.1 9)  leads  to 

°  i.N-1  “n.N-I  “  V  i.i  * 


°i  <  °  i,N-l  EN. 


where  r^  ■  Ms*  rj„1/2* 


.  -1  -1  . 
N-l  rM(“  i.i  “  “  i,M-r  * 


(3.27) 


Vl  Vl  ^Vl 

'**"1  2  Logn^j  (Logn^j)2 


end 


1  .  1,.  ,  ,  „  ,  Vl 

"  4  -  '  4(1"Vl)  *  *N,N-1  ~ 


3/4 


(1  -  ^  <  0. 


Hence,  (s  +  |)  1  1#  ^  ^  2  •  2£  n~1  beclltt*e  (3.18).  end  the 
inequality  (3.27)  iaplies 


-  (r^  +  ^)  £  (a  +  |)  cfj  <  tjj. 


Froa  (3.23),  0  i  j  j  if  0<Uj£I.  Finally,  we  obtain 


llA_1J(Q(u))l  -  Snp  lA~xJ(Q(u))u|  ± 


-1, 


Uni-1  -«(.+l/2>* 

Therefore,  the  following  two  Theorems  are  resulted: 

Theoren  3^:  The  siapping 


P(u)  ■  A-1  Q(u) 


upper  bound 

(3.28) 


(3.29) 


is  contract iwe  if 


a  - 


(3.30) 


where  the  natriz  A  and  the  vector  Q(u)  are  defined  in  (3. IS).  In 
another  word,  P(u)  has  unique  fized  point  in  this  case. 


Theorea  3.3:  When  h  <  2s.  the  senil inear  scheme  (3.10)  has  unique 


20 


solution,  and  it  e«n  be  solved  by  the  following  ' simple'  iteration 


0  (0)  - 


(3.31) 


d  +  Q  <  U(k_1)  )  .  k  >  1. 


Furthemore.  froa  (3.4).  U  is  a  aonotocally  increasing  sequence 


of  j  and  Q(u  )  is  nonnegative.  In  fact 


Q(pj+1-l,pj-l.pj_1-l)  -  Q(pJ+1,pj.pj_1)  +  Q* 


where  Q*  -  (pJ+1~Pj )  Cr£— -  “  - rrj - 7~  )  - 

^  Log[(pJ  1-l)/(pj-l)) 

(pj-pj_1)  - j  " - rrr  ) . 

Logp  Logl(pj-l)/(pj  -1)] 

Froa  Lena  3.3  (ii).  0(pj+1  ,pj ,pj_1)  >  0.  and  sinoe  p  >>  1 


Log[(p^+1-l)/(p^-l)]  -  Logp  -  p 


1/2  -1/2.2 

_ 


By  induction,  it  can  be  proved  that,  for  this  Model  problea  (3.1),  ^ 
in  the  iteration  (3.31)  preserve  the  aonotonically  increasing  with  j  for 
each  k.  Hence,  when  k  tends  ",  the  aonotonically  increasing  property 
still  reaain  true,  i.e.. 


CjiDj+r  j  “  O’1’2'*-- 

Now  we  consider  the  convergence  in  another  Banning:  the  oonvegence  of 
the  solution  of  the  nonlinear  difference  equations  (3.10)  to  the 
exact  solution  u  of  the  differential  equation  (3.1).  In  aatrix  fora 
(3.1)  with  the  local  truncation  error  foxaula  (3.12)  or  (3.13)  can  be 

expressed  as 


A  u  ■  d  +  Q(u)  +  Tr(u) 


(3.32) 
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Subtracting  (3. IS)  froa  (3.32)  yields  tbe  error  equation: 

A(n  -  U*)  «  Q(n)  -  (MU*1)  +  Tr(u).  (3.33) 

Vben  h  <  2e,  A  1  exists,  the  above  f omnia  is  equivalent  to 

u  -  U*  -  A_1(Q(u)  -  Q(Uh))  +  A_1Tr(u) .  (3.34) 

According  to  the  Lena  3.4,  for  h  <  2s,  A  *Q( u)  is  a  contraction  napping 
in  the  aaxiaua  non,  hence,  froa  (3.28) 

Ha-1(Q(u)  -  (MU*))!.  i  4(-"172)in  "  A.' 


Applying  Lena  3.1  and  (3.13),  we  obtained 

■»  -  0*1.  <  l  1  -  iTfTi75)  l'1  <  c  h2 

where  the  constant 


c  -  f-.-1/(2*> 


-3 

e 


3 


<  3.6  . 


(3.35) 


Therefore,  ve  have  proved  the  following  aain  result  of  this  section  : 


Thcorea  The  solution  of  the  seai-linear  scheae  (3.10)  converges 

to  the  exact  solution  of  the  singularly  perturbation  boundary  value 
probleu  (3.1).  The  rate  of  convergence  in  the  aaxiaua  non  is  of  second 
order.  Moreover, 

R  Uh  -  u  im  <  Ch2.  if  h  .<  2e .  (3.36) 

where  the  coefficient  C  defined  in  (3.3S)  is  uniforaly  bounded  for  all  s 
and  C  tends  to  zero  as  s  does. 


Note  that  the  restricted  aesh  condition  h  <  2s,  caused  by 


introducing  the  seai-linear  scheae,  is  only  used  in  the  steeper  gradient 


interval  z  >  1/2.  Hence,  it  is  possible  to  restrict  tbe  aesh  condition 
only  in  tbe  interval  .  As  a  natter  of  fact,  in  block  natriz  fora,  tbe 
acbeae  (3.10)  can  be  written  as 


I  A1t  A.,  I  IVl 

|  11  “12  i,  , 

•  *21  *22  1  101 


I  0  I  I  0  | 

I  I  +  I  I 

I  d  I  I  Q(U.V)  I 


Let  b^  and  b  be  two  different  uniform  aesh  sizes  used  in  tbe 
intervals  z  <  1/2  and  z2  1/2,  respectively,  and 


Applying  Gaussian  eliaination  of  aatriz  fora,  tbe  above  syatea  can  be 
rednoed  to 


A^O  -  d  +  (HD.-A^A^D) 


(3.37! 


where 

a  .  ,  *  x  •  .  .  1  *  *+1/2 

*22  (  *  ij  ’’  °  11  *  2  a^+572' 

o*^  -  Gy  if  (i,j)  *  (1.1). 

Henoe,  when  a  2  »j.  i.e.,  b  1  hj,  tbe  Lease  3.1  still  bolds,  tbe  above 
conclusion  can  be  iaproved  further  by  tbe  nezt  assertion. 

Xhmw  2,2:  If 


h  £  ain  (  bj,  2s  )  (3.38] 

then  the  iterative  procedure  (3.31)  converges  for  tbe  syatea  (3.37)  and 
tbe  error  estiaate  (3.36)  is  still  valid. 
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4.  Linear  second  order  two-point  boundary  layer  probleaa 

Consider  a  general  linear  second  order  singular  perturbation 
problen 


Lu  s  -e  u"  +  p(x)  u*  +  q(x)  u  “  f(x),  in  [0,1] 

u(0)  -  u(l)  -  0  (4.1) 

where  t  is  a  snail  positive  paraneter  and  p(x),q(x)  and  f(x)  are 
sufficiently  suooth  that  their  second  derivatives  are  uniformly  bounded 
for  all  x  in  [0,1]  and  for  all  e  >  0,  besides,  p(x)  2  p  >  0  ,  q(x)  2  0 
and  q(x)  -  p'(x)/2  2  6  >  0  on  [0,1].  Let 

LjU1^  -  -(o+Pj/aj^.j  +  tte+qjh2)!^  +  (Pj/2  -  sOU*^.  (4.2) 

The  interval  [0,1]  is  divided  into  tjro  subintervals  :  [0,1]  *  Ir  +  1^, 

where  1^  is  called  a  regular  subinterval  over  which  the  first  derivative 
of  u(x)  is  bounded  by  a  control  number  |»h,  and  I#  defines  a  singular 
subinterval  over  which  u'(x)  may  be  very  large. 

U  -U 

Define  p  ■ * —  n  t .  Similar  to  (ref[3.10]),  in  this  case  the 

i  2h 

semi-linear  scheme  becomes 

hfjh  if  Itl  £ 

LjU*  -  {  (4.3) 

h£*  -  iPy  ^j+l>  if  \x\  >  ph 

where 


■  r(.,c,b+c)  -  *§*  -  Ci,U*b'n*.>V 


(4.4) 


h 


I 


h 

\ 

l 

[ 

l 

h 


-  > 


4 


4 


< 


4 
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and  c  ic  a  parameter  to  be  chosen.  The  purpose  of  introducing  c  is 
two-fold.  First,  to  sake  the  scheae  to  be  well  defined;  second,  to  lead 
to  a  better  approxiaation. 

Now  the  corresponding  aatrix  A  in  (3. IS)  becoaes 


2a+hq2 

-<a+p2/2) 


-(a-Pj/2) 


2a+hq2  -(a-pj/2) 


~{a+pN-2/2)  2*+hq^_2  -(•-Pn_2/2) 

~<a+pN~l/2)  2a+hqN-l 


(4. 


Denote  the  detexainant  of  the  first  j  and  the  last  N-i  principal 


deterainant  by  D.  and  D,  „  , ,  respectively, 

j  l#  fr'l 

„  _  Dn-1  „  _  °j +1,^-1 

P-  - »  Pi  u_i  -  * 


then 


»  0  o 

n  "j ,N~1 


P__ ,  P_ 

-  (2a  +  hq  )D  -  (a  -  -2-±) (a  +  —  )D  ,, 
n  n  n-A  A  «  n*"A 

2  2 


DQ  *  1,  Dj  -  2a  +  hq2< 


K 


2a+hqn  -  (a-pa_1/2)(a+pa/2)pn_1 
Dj,N-l  “  <2a  +  hqj)Dj+l,N-l  U  "  ,(a  + 


°N,N-1  *  1#  DN-1,N-1  “  2a  +  hqN-l* 


(4.6) 


r  i  u.1  • 

j>  2..^  -  <*-Pj'2M**pJtl/2>f>J+liN_l 

Leaaa  4 .1 :  Assnae  that 


(i) 


(4.7) 


(ii>  qj  ^  2h<pj+l  "  Pj-1)'  for  aU 


(4.8) 
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then 

*,  *  ^72  •  f“  *u  *  <  x-1- 

s».n-i  *  I7T73  ■  £“  *u  *  <  "-1-  <«•»> 

fi 

Reaark:  At  a  discrete  fora  of  an  inequality  q(x)  2  p'(x),  (4.8)  it  a 
sufficient  condition  of  an  elliptic  fora  for  the  equation  (4.1). 

fissS ■ 

The  first  inequality  (4.9)  is  trivial  for  n*l.  By  induction, 
suppose  (4.9)  holds  for  n-1 ,  then  froa  (4.6)  with  the  condition  (4.8) 

*  lu  *  5(>Vi  "  Vx>  ‘  (*-Vi/llf‘  ■  “  * 

The  seoond  part  of  (4.9)  can  be  proved  in  the  saae  Banner. 

Definition:  f  real  n  x  n  aatrix  A  ■  (  o  .  )  with  o  .40  for  all  i#j  is 

ip  J  J 

called  a  G-aatrix  if  A  1  2  0  and 

•\j  2  11,21 

1  o’1,  ^  if  i  <  j  (4.10) 

Reaark:  A  is  a  G-aatrix  aeans  that  A  1  is  a  good  discrete  Green's 
function  in  the  aachanical  sense. 

Theorea  4.2:  The  aatrix  A  defined  in  (4.5)  is  a  G-aatrix  if  (4.7)  and 
(4.8)  hold. 

Pyoof : 

In  fact,  we  have 

i-1 
TT 
k-j 


i.j  '  i+l,N-l 


j-1  p"  D 

TT  (.  -  — )  D  — 

k«i  2  i  1,N^1  DVI  , 


if  i  <  j. 


Hence,  in  the  case  of  i  >  j 

e-1 .  .  D.  ,  1 

_ ii±_  .  - ....  Jr1 _ _  _ 

and  in  the  case  of  i  <  j 


i,J-l 


°  i.j  (*"Pj/2)Dj+l.N-l  (t+pj/2)Pj,N-l 


2  1.  O.E.D. 


Definition:  In  a  real  n-diaensional  apace  Rn  a  sapping  Q(x)  it  called 
semi-linear,  denoted  by  Q  s  SLn,  if  Q  e  C1  and  there  exists  a  constant  I 
snch  that 


{J<Q(x))(x+K))j  -  Qj(x) ,  for  1  <  j  <  n. 


(4.12) 


where  J(Q)  is  the  Jacobi  satrix  of  Q. 


In  particular,  any  linear  napping  Ax  is  semi-linear  for  every  square 
matrix  A  when  1*0.  Bence,  the  sen i- linear  sapping  defined  by  the  above 
definition  is  an  extension  of  the  linear  transformation. 


It  is  obvioas  that  SLn  forms  a  linear  space,  i.e,  if  both  Q(x)  and 
P(x)  are  semi-linear,  and  a.b  are  constants,  then  aP(x)+bQ(x)  is  also 
semi-linear.  Furthermore,  if  a  and  b  are  constant  vectors  with  a*(a^), 
b“(bj),  then  the  sapping  (SjPj(x)  +  bjQj(x))  is  also  semi-linear. 


As  an  example. 


R(x)  -  (RAx)),  RAx) 


Vi  '  xj 


Log((c+x. . , )/(c+x. ) ) 
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is  seni-linear,  because 


— *Lx  +  ■  ** — s . , ,  -  K. . 

a_  i  J+1  J 

dxj  dxj+l 

Moreover,  setting  s  shift  transformation 
matrix  form  the  system  (4.2)  becomes 


U*.  +  c  (1  i  j  i.N-1),  in 


A  V  -  d  +  Q(V) 


vbere  A  is  defined  in  (4.5),  (L  *  0  or  (r(Vj+1,Vj )  -  r(Vj-Vj_1) ) , 
depending  upon  whether  jelr  or  jel#,  respectively.  Since  R(x)  is 
semi-linear,  we  have 


{J(Q(V))V}  -  (Q(V)} if  j  <  N-l. 

J  J 


(4.13) 


Theorem  4.3:  For  given  g^(x)  defined  in  (4.4)  and  n^wn^wO,  the  mapping 


Q(x)  -  (gj(x))  is  semi-linear. 


Sinoe 


where 


•l  *  IA'lj,<>,V)IVI1  -  ■  “N.N-l  ' 

Nil  _! 

j-  n  {  Pjfl  i,j  "  Pj"lB  i'J“1  "  P»“  i.nrn-l/2 


Cj-l/2  *  r(Vj'Vj-l)*  *N,N-1  "  rN-l/2  + 


N-l/2 


it  is  easily  seen  that  if  the  condition  (3.30)  is  changed  to  (4.7)  ,  the 
Theorem  3.4  will  remain  true  when  p*20.  Therefore,  we  obtained  an 
extension  of  Theorems  3.4  and  3.5: 


Tb*orem  4.4:  Suppose 


p'(x)  2  0,  and  qj  2  '  for  *n 


(4.14) 
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l 

I 

! 

A 

1 

I 


A 


4 


I 


4 


then  the  napping  P(V)*A  *Q(V)  i$  contractive  if  (4.7)  holds. 

When  p'(x)  is  negative  sonewhere,  the  derivation  of  conditions  for  the 
napping  P(V)*A  1Q(V)  to  be  contractive  is  little  nore  conplicated.  In 
this  case,  set 


where  the  two  terns  on  the  righthand  side  denotes  a  positive  parts  and  a 
non-positive  part,  respectively.  Suppose  (KV^il ,  since  O^rjlr^,  fr on 
the  Theoren  4.2 


,i  '  Pj°  l.j  *  Pj-l“  i.j-l  *  V  i.»Vl/2 

J  n 

1  ♦  v'uVi/i  * 

Ipj-1°  i.j-l  '  pj*  1. j1  lr«“rj-l/21 

-1  -  + 

<  v  i.i‘  Pi  *  ^  (pj-rpj)  >• 


Hence 


-»i  <  !  *  2(.ipl+1/2)  /J  Vl-'j1' 

In  the  sane  nanner  of  the  derivation  in  the  last  section,  we  obtain  the 
following  result. 


Theoren  The  conclution  in  Theoren  4.4  still  holds  even  if  the 

condition  p'2,0  is  renoved  provided  that  the  iequality  (4.7)  is  changed 
to 

Bp*.  i 

a  >  Mart  - ,  2  1  (p..-p. >  -  SA,! •  (4.15) 

2  j-n  j  1  j  2  N 

In  order  to  get  an  error  estinatation,  using  (3.34),  it  is  only 
need  to  find  a  bound  for  |A_1Tr(  u)Hs  when  the  napping  is  contractive. 
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This  reduces  the  error  estimates  to  estimate  the  bond  of  the  trncatioa 
error.  Substituting  the  exact  eolation  a  in  the  scheme  (4.3)  yields 

Vj  ■  “j  * 

if  j  e  Ir*  and 

Li»j  -  “j  '  •j'Vr'j'Vi’  *  V’' 

if  j  e  I#,  where  0<{,q<l. 

It  has  been  noted  that  the  exact  aolntion  has  a  factorization  which 
consists  of  two  parts  :  one  is  regular,  another  is  singular  (see  Hon 
Honde[3],  Lemma  1),  hence 

u(x )  -  T{*(x)  +  Z(x)} 

T(x)  -  e~p(1)(1"x)/*  -  x  -  (l-x)e"p(1)/#  (4.17) 

where  y  ■  1  in  1  in  euf(x)/p(l)  is  a  constant  bonded  niformly  for  all 

a— >0  x— >1 

0<e<l,  and 

|Z(x)  I  i  C.  |Z*  (x)  I  i  C,  |Z' '  (x)  I  ±  C{l+7  (4.18) 

e 

C  is  a  constant  independent  of  e,  and  0<P.£p  . 

Set  c  in  (4.4)  eqnl  to  y  which  can  be  also  obtained  in  computing 
process.  Because  the  singularity  of  the  exact  solution  u(x)  is  only 
near  x*l,  the  width  of  the  bondary  layer  in  which  u(x)  has  large 
derivatives  is  less  than  k  times  s,  where  k  is  a  constant  no  natter  how 
small  s  is,  this  implies  that  we  can  take  I#  ■  (l-ke,l),  and  on  [0,1-ke] 
u(x)  and  its  first  fourth  derivatives  are  niformly  bonded.  Hence  we 
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only  seed  to  consider  the  aaxiaua  error  within  the  boundary  layer.  Froa 

(4.17).  when  h  £  eflpfls/2,  in  the  interval  I# 

3 

(4)  (3)  .  p(1)  .  ,  .  M x,  -p(l) (l-x)/s 

-eu  +  pn  *  c{  —  j—  lp(x)-p(l) Je 

eZ(4)  +  pZ(3)  )  -  0(h"2), 


...  -iC«- Lf... 


C+U  B+C 
K  ^ i. vl 


(u* ' (n+c)  -  n'2) 


.  <»-*>/• .  z..  a-*)/. .  2  ,  cJ 


n+c  “  2 


-  tI>(l)e-p(l)(l-x)/.  +  z]2} 


.  -  2]  ♦  ... 

n+c  •  • 

(n"  -  sd}'  «  £L,“*,<1)(1”*)/,(Eili)2[Eii>(z+l)  -  2)  ♦  ... 
c+n  n+c  e  c 


Since  n(l)-0.  Z(l)-  -1,  eo 


i"-^-°<h~1).  (u" 


n*2  -2 

kJ*  “0(h  > 


(4.19) 


Substituting  into  (4.16)  yields 


aaxlTr(x) I  ■  0(h) . 

I 

s 

However,  for  the  nsnal  linear  central  difference  scheae  the 
corresponding  truncation  error  is  O(h^).  Hence,  introducing  the 
seni-linear  scheae  gives  one  aore  order  of  precision  for  the  aaxiana 
local  truncation  error.  Considering  that  the  width  of  the  boundary 
layer  I(  is  only  he,  the  nuaber  of  knots  using  the  seai-linear  scheae  is 
alwsys  less  than  a  constant  if  the  ratio  ^  reaains  a  constant.  Bence 


lA^TrU)!,  -  0(h) 


(4.20) 


Since  the  error  systea 


31 


A  (u  -  U*)  -  Q(u)  -  Q(U)  +  Tr(u) , 

so 

No  -  U^l  i  lA~1(Q(*)-Q(Uk»l  +  BA_1Tr(x)||, 

No  -  U*N  <  {  1  -  |A_1J(0(o))N  }_1Ia""^t*(x)N. 
if  Ba_1J(Q(o))H  No  -  u^l  <  1. 

Therefore,  we  hive  extended  the  error  estimate  Theorea  3.6  for  the 
general  problem  (4.1). 

Theoren  4.6:  The  solotion  of  the  semi-linear  scheme  (4.3)  converges  to 
the  exact  solotion  of  the  singolar  perturbation  problem  (4.1)  as  h  tends 
zero  if  the  (4.8)  holds.  Moreover,  if  the  ratio  of  h  to  s  is  less  than 
a  constant,  as  given  in  (4.7)  or  (4. IS),  aeeording  to  the  condition  of 
Theorems  4.4  and  4.S  respectively,  then  there  exist  an  error  estimation 
such  that 

N  U*  -  o  I.  <  Ch.  (4.21) 

where  the  coefficient  C  is  uniformly  bounded  for  all  e>0. 

Using  the  same  block 'matrix  technique  described  in  the  last  section,  the 
mesh  constrain  (4.7)  or  (4. IS)  can  be  limited  only  in  the  interval  1(. 

In  the  regular  interval,  a  larger  mesh  step  can  be  allowed. 


t _ _ _ 


5.  Computational  r« salts 


In  this  section,  three  numerical  examples  are  presented  to  show  the 
effect  of  the  semi-linear  scheme.  Only  uniform  mesh  sizes  are  need 
thronghont  this  section.  The  emphasis  is  on  the  comparison  among  the 
semi-linear  scheme,  the  corresponding  linear  centered  difference  scheme, 
and  an  'Upwind  method'  described  by  Christie  and  el.  in  [3]. 

EXAMPLE  1.  —  Model  Problem 

Ln  ■  -an"  +  u'  ■  0  in  (  0,1  ) 

u(0)  -  0.  u(l)  -  1  (3.1) 

Nnmerical  results  for  the  model  problem  (3.1)  are  tablated  in  Table  1-3, 
and  a  comparison  of  throe  different  methods  is  given  for  s  “  1/60  and  h 
■  1/40  in  the  sense  of  the  pointwise  error.  For  the  linear  centered 
difference  scheme,  there  is  a  relative  error  of  at  the  last  Interior 
mesh  point  z  “  0.975.  The  'Upwind  Symmetric  Quadratics  Method'  has  a 
relative  error  1#2%  at  the  same  point.  The  advantage  of  the  semi-linear 
scheme  is  clear.  Three  or  seven  digits  can  be  improved  for  an  iteration 
error  of  0.1D-3  or  0.1D-8,  respectively.  For  a  given  accuracy,  the 
required  CPU  time  nsing  the  semi-linear  scheme  is  twice  less  than  using 
the  usual  linear  scheme. 

For  the  same  s.  Table  2  describes  the  convergence  rate.  If  the 
given  admissible  error  is  0.SD-4,  the  semi-linear  scheme  requires  a  mesh 
size  h  approximately  1/30,  while  the  linear  scheme  requires  a  mesh  size 
1/2000.  The  total  CPU  time  between  these  two  schemes  differed  by  a 
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factor  of  2  to  3.  If  we  consider  tbs  error  of  the  first  derivative,  the 
difference  between  the  two  scheaes  is  even  greater  because  we  nse  a 
sen i- linear  differentiation  foraula  to  aatch  the  seai-liaear  scheae 
(3.10)  in  order  to  get  the  discrete  first  steep  derivative. 

The  influence  of  the  the  convergence  of  the  aeai-linear 
discrization  is  presented  in  Table  3,  depending  upon  the  ratio  of  t  to 
h.  Ve  consider  three  cases:  h/e  “  2,  3/2,  and  1.  According  to  Theorea 
3.2,  the  iteration  converges  if  h/e  £  2.  These  nvaerical  results 
satisfy  the  theory.  The  iteration  is  also  convergent  for  h  ■  2s, 
however,  this  is  not  recoaaended  because  the  rate  of  convergence  is  too 
slow.  For  a  practical  choice  of  h/e  ■  1.3,  a  higher  precision  seeas 
with  obtainable  less  CPU  tiae. 

Hence,  to  get  the  saae  accuracy,  one  can  solve  a  saall  seai-linear 
systeas  instead  of  the  original  large  linear  systea.  However,  the  error 
in  the  two  scheaes  behaves  differently.  For  the  linear  scheae,  when  the 
ratio  a/h  is  constant,  the  error  is  also  constant  independent  of  the 
size  of  e.  For  the  seai-linear  scheae  the  error  converges  to  zero  as  s 
does.  The  ezaaples  to  follow  will  show  the  saae  phenoaienon. 
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Table  1. 

(  e  «  1/60,  h  ■  1/40  ) 


X  ( 

Theoretical 

Linear 

Semi-linear 

Upwind 

! 

solution 

scheme 

ei  -  0.1-3 

ei  ■  0.1-8 

method 

0.800  j 

0.000006144 

0.000000173 

0.000006141 

0.000006144 

0.850  | 

0.000123410 

0.000008500 

0.000123447 

0.000123410 

- 

0.875  I 

0.000553084 

0.000059499 

0.000SS3305 

0.000553085 

- 

0.900  | 

0.002478752 

0.000416493 

0.002477749 

0.002478752 

0.0026 

0.925  I 

0.011108997 

0.002915452 

0.011110675 

0.011108997 

0.0115 

0.950  | 

0.049787068 

0.020408163 

0.049785741 

0.049787068 

0.0510 

0.975  I 

0.223130160 

0.142857143 

0.223130576 

0.223130160 

0.2258 

1.000  1 

1.0 

1.0 

1.0 

1.0 

1.0 

Max  error  of  n 

-0.8027-01 

0.1679-5 

0.2044-9 

0.27-2 

Location 

of  max  error 

0.975 

0.925 

0.875 

0.975 

Max  error  of  u' 

0.56+1 

0.95-4 

0.11-7 

- 

Number  of  iterations 

1 

12 

22 

1 

CPU  TIME 

(seconds) 

0.07 

0.15 

0.26 

— 

Remark:  ei  —  the  admissible  error  for 

iterations  of 

the  semi-lineat 

scheme . 

Table  2 


(  e  *  1/60, 

ei(SLl)  - 

0.1-4, 

ei(SL2)  -  0.1 

-8  ) 

1/h 

M 

20 

30 

40 

60 

640 

1920 

Max 

L 

-0.25+0 

-0.14+0  -0.80-1 

-0.35-1 

-0.27-3 

-0.30-4 

error 

SL1 

-0.46+1 

0.54-5 

0.17-5 

-0.35-6 

-0.37-9 

-0.51-12 

of  u 

SL2 

-0.46+1 

0.32-9 

0.20-9 

-0.51-10 

0.46-12 

0.46-16 

Max 

L 

0.70+1 

0.66+1 

0.56+1 

0.39+1 

0.80-1 

-0.25-1 

error 

-0.10+4 

0.31-3 

0.95-4 

-0.21-4 

-0.13-6 

0.50-9 

of  u' 

SL2 

-0.10+4 

0.19-7 

0.11-7 

-0.28-8 

0.16-9 

0.51-13 

CPU 

L 

0.01 

0.02 

0.07 

0.08 

0.46 

1.30 

SL1 

0.93  (100) 

0.51  (35) 

0.15 

(12)  0.16  (7) 

0.90  (3) 

2.78  (3) 

TIME 

SL2 

0.93  (100) 

0.74  (71) 

0.28 

(22)  0.24  (12) 

1.10  (4) 

3.21  (4) 

Remark.  1.  The  number  in  the  bracket  is  the  number  of  iterations  required  to 
reduce  the  error  to  be  less  than  the  admissible  range  ei. 

2.  Max  o'  “  60. 


« 


1 


I 


Table  3-1 


I  I  l/s  1 

i  10 

20  < 

1  ei 

1  0.1-4 

0.1-4  0 

Error 
of  n 


Error 
of  o’ 


L  I  < -  -0.14+0  - > 

SL  I  -0.45+0  -0.95-4  0.32-9  0.20-12 


0.01 

0.21 

(100) 


-0.20-2  0.19-7  0.18-10 


3 


Remark.  1.  The  number  in  the  bracket  ia  the  number  of  iterations  required  to 
reduce  the  error  to  be  less  than  the  admissible  range  ei. 


Table  3-2 


h/e  -  3/2 


Error 
of  n 


I  < - 

I  -0.23-3 


Error  1 
of  n*  1 

1  L 

1  SL 

1  0.14+1 

I  -0.39-2 

CPU  1 

1  L 

1  0.01 

| 

TIME  j 

1  SL 

1  0.03 

f 

1 

1  (  12  ) 

—  -0.80-1  - 

0.15-5  0.20-9 


.28+1 

.43-4 


0.1-12 


- > 

0.28-13 


0.24-11 


Table  3-3 


r  i 

l/e  1 

1 

ei  1  0.: 

Error 
of  n 


Error 

1 

L 

of  u' 

1 

SL 

|  < - 

I  -0.49-2 


I  0.01 

I  0.01 

I  (  7  ) 


—  -0.35-1  - 

-0.92-3  -0.64-5 


0.01 
0.02 
(  7  ) 


-0.15-3 


0.02 
0.05 
(  7  ) 


-0.51-10 


-0.28-8 


0.1-12 


- > 

-0.69-14 


-0.59-12 


Figure  2 

THREE  METH0DS  C0MPARING  (EPS=l/60.  H=l/40) 


0.906 

0.938 

X 

1 

SOLUTION  X 

=  0.975*  U 

=  0.2231 

2 

LINEAR  SCHEME  U  =  0. 

H29 

3 

SEMIL INEAR  SCHEME  U  = 

0.2231 

4 

UPWIND  U  = 

0.2258 

Figure  3 


TABLE  3-2,  H=1.5EPS,  1/EPS=15. 30,60, 90 


1.0E+02  r 
l.QE+Ol  k 
1.0E+00  t 
1.0E-01  §- 
1.0E-02  j- 
1 . 0E-03  k 
1.QE-0*  S’ 
1.0E-05  6- 
l.QE-06  k 
1.0E-07  f- 
1.0E-08  \ 
1.0E-09  §• 

i.oe-io  r 


E- 


1.0E-12  k 

l.QE-13  S' 

l.OE-H  t- 

10.0 


32.5 


55.0 

X  =  1/H 


77.5 


LINEAR  SCHEME  FUNCTION  ERROR  =  0.081 
SEMI LINEAR  SCHEME  FUNCTION  ERROR  0.23E-3 
LINEAR  SCHEME  FIRST  DERIVATIVE  ERROR 
SEMI LINEAR  SCHEME  FIRST  DERIVATIVE  ERROR 
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H 


■ 


I 


I 


1 


« 


EXAMPLE  2.  A  linear  singular  perturbation  problem  with  constant 
ooeff icients 


Lu  *  -  an"  +  o'  +  (l+e)u  *  f(x),  in  (  0,3  ) 


n(0)  *  n(l)  *  0  (4.1) 

where  f(x)  *  (l+e)(a-b)x  -  sa  -  b,  a  ■  1  +  e  ^+e^e,  b  *  1  +  e  *,  with 
exact  solution 


,  .  -(1+e) (l-x)/e  -x 

o(x)  *  e  +  e 


a  +  (a-b)x 


It  is  not  difficulty  to  find  that  the  coefficient  y  in  (4.17)  here  is 
equal  to  1  without  knowing  the  exact  solution  in  advance.  Hence,  we 
assign  the  constant  c  ■  1  using  the  semi-linear  scheme  (4.2). 


The  result  listed  in  Table  4  shows  that  iterations  converge 
monotonically  if  the  ratio  h It  £  2  and  that  the  calculated  results  agree 
well  with  the  theoritical  analysis.  For  the  same  mesh  size,  required 
CPU  time  using  semi-linear  scheme  is  little  more  than  using  the  usual 
linear  scheme  (  about  10%  to  20%  for  small  e  ),  but  more  than  one 
significant  digit  is  obtained.  More  small  t  is  contained,  more 
advantage  the  semi-linear  scheme  has.  Therefore,  with  the  same 
accuracy,  using  semi-linear  scheme,  an  large  linear  system  arisen  from 
the  linear  scheme  is  replaced  by  a  smaller  semi-linear  system.  For 
instance,  when  e*0.01,  the  maximum  error  using  the  semil inear  scheme 
with  N*60  (after  13  iterations)  is  less  50%  than  of  the  usual  linear 
scheme  with  !M00,  meanwhile  the  ratio  of  CPC  time  is  0.27  :  1.12(sec.). 
Vhen  e*0.001  ,  the  maximum  error  using  the  semilinear  scheme  with  N*600 
is  less  ten  times  than  of  the  usual  linear  scheme  with  N*2000,  and  the 
ratio  of  CPU  time  is  3.22  :  7.87(sec.). 


Table  4-1 


e  “  0.1  ei  ■  1.0-5 


1 


3  9 


[• 


h-l/N 

z 

Max(Er(u) ) 

Max(Er(n'>) 

CPC 

NE 

5 

L 

0.800 

0.1935D+00 

0.124704-01 

0.01 

1 

SL 

0.800 

0.1935D+00 

0.1247D4-01 

0.01 

2 

6 

L 

0.833 

0.1474D+-00 

0.116  304-01 

0.01 

1 

SL 

0.833 

0.4765D-01 

0.8202D-01 

0.03 

10 

10 

L 

0.900 

0.6548D-01 

0.7837D4-00 

0.01 

1 

SL 

0.800 

0.2255D-01 

-0.264204-00 

0.08 

7 

20 

L 

0.950 

0.1974D-01 

0.3247D+00 

0.02 

1 

SL 

0.850 

0.6793D-02 

-0.111204-00 

0.04 

4 

40 

L 

0.975 

0.5493D-02 

0.1057D4-00 

0.05 

1 

SL 

0.850 

0.1499D-02 

-0.3305D-01 

0.10 

4 

Remark. 

1.  NE 

is  the  amber 

of  iterations  required  to 

redace  the  error 

to  be  less  than 

the  admissible  range  ei. 

Table  4-2 

c  -  0.01, 

ei  ■  1.0-5 

h»l/N  | 

X 

Nax(Er(u)) 

Max(ErU')) 

CPC 

NE 

50  I 

L 

0.980 

0.1595D+00 

0.111604-02 

0.15 

1 

SL 

0.960 

0.795 8D-02 

-0.182904-01 

0.40 

24 

60  | 

L 

0.983 

0.123 OD+OO 

0.1020D4-02 

0.18 

1 

SL 

0.967 

0.6016D-02 

-0.154304-01 

0.27 

13 

100  I 

L 

0.990 

0.55610-01 

0.658104-01 

0.29 

1 

SL 

0.970 

0.28720-02 

-0.618404-00 

0.43 

7 

200  ! 

L 

0.995 

0.16880-01 

0.262404-01 

0.57 

1 

SL 

0.970 

0.63390-03 

-0.174204-00 

0.69 

4 

400  I 

L 

0.998 

0.47050-02 

0.836404-00 

1.12 

1 

SL 

0.970 

0.16810-03 

-0.57430-01 

1.36 

4 

i 


40 


Table  4-3 


e  »  0.001, 


•i  *  1.0-5 


b-l/N  I 


Naz(Er(n) ) 


Max(Er(o' ) ) 


0.998 

0.992 

0.998 

0.995 

0.999 

0.995 

0.999 

0.995 

0.999 

0.996 


0.15680+00 

-0.1981D-02 

0.1210D+00 

0.93130-03 

0.5475D-01 

0.2977D-03 

0.4045D-01 

0.21500-03 
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Figure  5 


FR0M  TRBLE  4  H/EPS=5/3  EPS=0. 1,0.01,0.001 


1  LINEAR  SCHEME  FUNCTION  ERROR 

2  SEMI LINEAR  SCHEME  FUNCTION  ERROR 

3  LINEAR  SCHEME  FIRST  DERIVATIVE  ERROR 

4  SEMILINEAR  SCHEME  FIRST  DERIVATIVE  ERROR 
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1 
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Figure  6 


FR0M  TABLE  4  CPU  C0MPARING 


PRECISI0N 


1  LINEfiR  CENTERED  DIFFERENCE  SCHEME  EPS=0.1 

2  SEMILINEAR  DIFFERENCE  SHEME  EPS=0.1 

3  LINEAR  EPS=0.01 

*  SEMILINEAR  EPS=0.01 

5  LINEAR  EPS=0.001 

6  SEMILINEAR  EPS=0.001 
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EXAMPLE  3.  A  seal-linear  singular  perturbation  problen 

Ln  ■  -  su"  +  n*  +  (l+e)u  ■  f(x,u),  in  (  0,1  ) 

u(0)  “  n(l)  “  0  (5.1) 

with  the  saae  solution  as  Exsaple  2,  where 

_  c 

f(x,u)  »  a  -  b  -  (l+e){e  x  -  u  +  .  ■  -  ■■  ) , 

u+s-(s-b)x-e 

s-l+e-<1+e)/‘.  b-l+e^.e  -  e2(1+*)(1“x)/‘. 

This  exsaple  shows  the  advantage  of  using  SL-  scheae  over  L-scheae 
wore  than  exsaple  2,  because  both  seheaes  have  to  be  solved  by  iteration 
for  this  seai-linear  problea,  and  in  this  ease  SL-scheae  requires  less 
CPU  tiae  than  in  L-scheae  for  the  saae  aesh  sire. 

Two  different  results  of  two  seai-linear  seheaes  are  listed  in 
Table  5.  SL1  seheae  is  described  in  section  4.  SL2  scheae  cones  froa  a 
Seai-linear  Gslerkin  aethod  [Jiaehang  Sun,  11].  Usin^  the  saae  siaple 
iterative  procedure  for  the  saae  size  h,  SL2  secas  need  less  auaber  of 
iterations  required  to  reduce  the  error  to  be  less  than  the  saae 
adaissible  range  than  SL1,  but  SL2  costs  s  little  aore  CPU  tiae  to 
coapute  its  coefficients.  At  any  rate,  the  advantage  of  both 
seai-linear  seheaes  over  the  linear  seheae  is  elesr.  The  stronger  the 
singularity  of  the  solution  of  the  differential  equation,  the  aore  the 
advantage. 
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14 
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0.25610-03 
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Table  5-2 

•i  -  1.0-5 

b-l/N 

X 

Max(Er(n)) 

UaxtErU')) 

CPU 

NE 

50 

L 

0.980 

0.13 91D+00 

0.11160+02 

4.33 

100 

SL1 

0.960 

0.8372D-02 

-0.18860+01 

4.72 

100 

SL2 
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Table  5-3 

e  -  0.001. 

ei  -  1.0-5 

h«l/N 
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Figure  7 


FR0M  TABLE  5  H/EPS=5/3  SEMI LI NEAR  EQUATI0N 


1  LINEAR  SCHEME  FUNCTI0N  ERR0R 

2  SEMIL INEAR  SCHEME  FUNCTI0N  ERR0R 

3  LINEAR  SCHEME  FIRST  DERIVATIVE  ERR0R 

4  SEMI LINEAR  SCHEME  FIRST  DERIVATIVE  ERR0R 


FUNCTI0N  ERR0R 


Figure  8 


TABLE  5-3  EPS=0.001  SEMI LI NEAR  EQUATI0N 


1  LINEAR  SCHEME  FUNCTION  ERR0R 

2  SEMILINEAR  DIFFERENCE  SCHEME 

3  SEMILINEAR  GALERKIN  METH0D 


CPU  SEC0NDS 


Figure  9 


FR0M  TfiBLE  5  CPU  C0MPARING 


PRECISI0N 


1  LINEAR  CENTERED  DIFFERENCE  SCHEME  EPS= 

2  SEMI LINEAR  DIFFERENCE  SHEME  EPS=0.1 

3  LINEAR  EPS=0.01 

4  SEMI LINEAR  EPS=0.01 

5  LINEAR  EPS=0.001 

6  SEMIL INEAR  EPS=0.001 
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